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Abstract 

A simple model of electron-vibron interactions in buckminster- 
fullerene ions is solved semiclassically. Electronic degeneracies of C6o n_ 
induce dynamical Jahn- Teller distortions, which are unimodal for n^3 
and bimodal for n = 3. The quantization of motion along the Jahn- 
Teller manifold leads to a symmetric-top rotator Hamiltonian. I find 
Molecular Aharonov-Bohm effects where electronic Berry phases de- 
termine the vibrational spectra, zero point fluctuations, and electrons' 
pair binding energies. The latter are relevant to superconductivity in 
alkali- fuller enes . 

PACS: 33.10.Lb,71.38.+i,74.20.-z 

For polyatomic molecules, the adiabatic approximation is often used to 
eliminate fast electrons in favor of an effective potential for the slow nuclei. 
This approximation requires special care when the positions of the nuclear 
coordinates are near points of electronic degeneracy. If the electron-ion inter- 
action is linear in the ionic displacements (a generic case for symmetric, non 
colinear molecules the classical Jahn Teller theory 0] predicts that the 
molecule distorts and some (or all) of the electronic degeneracy is lifted. The 
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classical JT theory is controlled by the largeness of S= \ E JT \/(huj), where uj 
is the characteristic vibrational frequency, and E JT is the relaxed energy of 
the distortion. For S = oo, (i.e. "strong coupling" or "classical" limit), the 
zero point motion of the ions is ignored. 

At finite S however, quantum corrections can be quantitatively and qual- 
itatively important. For example, while the JT distortion may break the 
Hamiltonian symmetry, quantum fluctuations along the degenerate manifold 
or tunneling between JT minima can restore the ground state symmetry 0. 
This is often called the "Dynamical Jahn Teller" effect. In addition, ion co- 
ordinates may be subject to quantum interference effects. Longuet-Higgins 
has found that a vibrational orbit which surrounds a point of two-fold elec- 
tronic degeneracy, can acquire a negative sign from the transport of the 
electron's wave function [[|. This effect, often dubbed as the "Molecular 
Aharonov Bohm (MAB) effect", produces in triangular molecules half-odd 
integer quantum numbers. This has been recently confirmed spectroscop- 
ically in Na3 0. The MAB effect is a simple example of the geometrical 
Berry phase, which appears in a wide host of quantum phenomena [EN. 

The soccer-ball shaped molecule Cqq (buckminsterfullerene) and its var- 
ious crystalline compounds have ignited enormous interest in the chemistry 
and physics community in past two years |§. Since the discovery of super- 
conductivity in A3C60 (A=K,Cs,Rb), with relatively high transition temper- 
atures (T c « 20°-30°K), much attention has been given to the electronic 
properties of charged Cgo™ - ions. Cgo a highly symmetrical molecule (a trun- 
cated icosahedron) , and its electronic lowest unoccupied molecular orbitals 
(LUMO) are three- fold degenerate. Thus the Ceo™ - ions are natural candi- 
dates for manifestations of dynamical JT effects and MAB effects discussed 
above. Several groups have identified the five-fold degenerate H g (d-wave 
like) vibrational modes that couple strongly to the LUMO orbit als0, [8|, §]. 
Varma, Zaanen and Raghavachari (VZR) proposed that these modes un- 
dergo a dynamical JT distortion and calculated the JT induced pair binding 
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energies at several fillings. These results were used to explain the large T c 's 
of fullerenes relative to doped graphite superconductors. VZR used the clas- 
sical approximation, and restricted their calculation to unimodal distortions 
(defined later). 

Density functional and deformation potential calculations for C6o _ @) !], 
estimate E JT « 40 meV. The important vibrations are in the range of 
hut ^0. 1-0.2 eV. Thus, the classical parameter is in fact quite small: S ~ 

0. 2-0.4. This indicates that the ions' quantum fluctuations cannot be justi- 
fiably neglected. 

In this paper the vibrations about dynamically distorted buckminsterful- 
lerine ions are quantized semiclassically. I extend previous work of O'Brien, 
who has solved the n = 1 case both exactly and semiclassically [0, |TTJ. First, 
the unrestricted classical JT distortions are determined. For C6o n ~, n^3, 
the JT distortions are unimodal, i.e. involve one quadrupolar mode in the 
principle axes frame. For Cgo 3- , the JT distortion is found to be bimodal, 

1. e. two modes are distorted simultaneously. Subsequently, the quantum dy- 
namics parallel and perpendicular to the JT manifold are determined. The 
excitation spectra and pair binding energies for n = 1, . . . 5 are determined 
up to second order in S 1-1 . I will show that that Berry phases give rise to se- 
lection rules for the pseudo-rotational quantum numbers. These kinematical 
restrictions effect the pairing interaction between electrons, and, therefore, 
also the superconducting transition temperature. 

This discussion is restricted to the simplest electron-vibron interaction 
model of Ceo, which captures the symmetries and degeneracies of this system. 
Electron-electron interactions are presently ignored. The wave functions of 
the LUMO t\ u states are represented by the L=l triplet \x), \y), \z). A single 
vibronic H„ multiplet is represented by five real coefficients 1T0, |7| , 
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where m, fi = —2,-1,... 2, and a\ m are the coefficients of the spherical 
harmonics Yi m ||12|| . The Hamiltonian is H = H el + iY 1 " 6 , where 
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where g is the dimensionless electron phonon coupling constant. H is 
invariant under rotations of q in R 5 and by construction, the eigenvalues of 
H are invariant under 0(3) rotations of the molecule's reference frame. H el 
is diagonalized by 
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where the Euler angles w = (0, 9, ip) define the 0(3) rotation which diago- 
nalizes H el . In the diagonal basis, the only non-zero vibrational components 
which couple to the electrons are q (°> = (r, 0, z, 0, 0). 

Since \q\ 2 is invariant under 0(3) rotations, the total adiabatic potential 

is 



Too 



V(z,r) ='^-( z 2 + r 2 ) + (ni(z- V3r) + n 2 (z + V3r) - n 3 2z) . (4) 



rii are the occupations of the orbitals = 1,2,3 (ordered from top to 
bottom in H el ), and J2i n i = n. V is minimized by the JT configurations 
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(z,r,fii), which yield the classical energies = V(z n ,f n ). The distortions 
are shown in Table I. We choose S = \g 2 as our semiclassical parameter. 

By ([I]), if we define axis 3 to be at 9 = 0, the z mode is described by 
z^(3cos 2 9 = 1), and the r mode by sin 2 9cos(24>). Thus by Table I, 
n — 1, 2, 4, 5 have unimodal distortions which are symmetric about the 3 axis, 
and n = 3 has a bimodal distortion, about the 3 and 1 axes. 

In order to quantize the vibrations, it is useful to express the kinetic 
energy in terms of small fluctuations about the JT distortion. To that end, 
we parametrize the JT degenerate manifold, {q^}, in terms of the Euler angles 
roofEq. (D: 

Q„M = E M^ m D%] m ,{w)M-^,qf . (5) 

m,m' ,fi'=— 2 

is the irreducible rotational matrix of angular momentum L |I2 |. The 
classical kinetic energy can be derived from (|5|) by the chain rule for differ- 
entiation. After some cumbersome, but straightforward, algebra we obtain 
the compact expression: 



1=1,3 

U\ = — sin ip6 + cos 9 sin 
u>2 = cos ipO + sin ip sin 9(p, 

C1J3 = ip + COS 9(f) 
(h,h,h) = ((^>z + r)\{VZz-r) 2 ,Ar 2 ). (6) 

For finite JT distortions, we can identify Ii(z, f) as moments of inertia in the 



principle axes frame fl2|| . Thus, the Euler angles dynamics follow those of a 



rigid rotator [13 



The unimodal and bimodal cases will be discussed separately. For the 
unimodal cases, f = and Ii = ^ 2 (3,3,0) on the JT manifold. The coordi- 
nate ip decouples from the rotational kinetic, which becomes that of a point 
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particle on the sphere. The quantization of the rotational part is therefore 
simply (2ii) -1 L 2 [|13|]. The remaining coordinates are three harmonic oscilla- 
tors r 7 = (r cos(2-0) , r sin(2-?/>) , z — z). Including the quadratic potential terms 
in V(z,r), we arrive at the vibrational eigenenergies: 

E™[L, n 7 ] = hu [^- 2 L{L + 1) + |> 7 + ±) j (7) 

The rotational parts of the eigenfunctions are 

n2(q) = Y Lm (6A) IIW ( g ) 

is 

where |rij S )' is an electron Fock state in the principle axes basis. The overlap 
of this Fock state with a Fock state of the stationary basis is a determi- 
nant that contains n factors of Y\ Vis - Under reflection, Yhmi^ — 0,<f) + n) ^ 
(— l) L YL m (8, 0). Therefore the electronic wave function yields a Berry phase 
factor of (— l) n for rotations between inverted points on the sphere which 
corespond to closed orbits of q. Due to the invariance under reflection of q 
(and thus the left hand side of (§)), a selection rule is obtained: (— l) L+n = 1. 
Thus, the ground state for n—1, 5 has pseudo-angular momentum L = 1 that 
contributes to the zero point energy. 

The analysis of the bimodal case of n = 3 proceeds along similar lines. 
From Eq.(|6|) and Table I, we see that (Ii,l2,h) = 3(? 2 (4, 1, 1). Thus, the 
kinetic energy includes the rotation of a rigid body with two equal moments 
of inertia. The quantization of this system is the quantum symmetric top 
Hamiltonian. Fortunately, its solution is a well-known textbook problem (see 
e.g. Ref.jjl], In addition to the rotator, there are two harmonic oscillators 
r 7 = (z — z, r — f). The eigenvalues of the bimodal C 60 3 ~ molecule are thus 
given by 

E U = ?*> U^L{L + 1) - ^ + |>7 + ~)) (9) 
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where L and k are quantum numbers of \L\ 2 and L\ respectively, and k < L. 
The rotational part of the eigenfunctions are 

*L*H = -dShIIW ( 10 ) 

is 

m is the eigenvalue of L z , where z is a stationary axis. In distinction to the 
unimodal case, there is no single reflection which fully classifies the symmetry 
of the wavefunction. However, one can obtain negative signs by transporting 
the electronic ground state in certain orbits. We define the rotation of tt about 
principle axis Lj as Cj. The Berry phases associated with these rotations can 
be read directly from Eq. (§). For example: for ip — > ip + tt (C 3 ), the states 
|1) and |2) get multiplied by (—1). Since D^' k transform as Y Lk under Cj, it 
is easy to determine their sign factors. The results are given below: 

C x : |1,0,2)' 
C 2 : 11,0,2)' - 
C 3 : |1,0,2>'- 

Clearly, q, being coefficients of quadrupole distortions, is invariant under 
Ci,C 2 ,C 3 . Thus, Ci describe closed orbits in R 5 . In order to satisfy ([[(]) and 
using the degeneracy of E bl for k —k, we find that L must be odd and k 
must be even. In particular, the ground state of (|9|) is given by L — l, and 
k = 0. 

A relevant quantity for superconductivity is the "pair binding" energy [0, 

HI 

U n = E n+ i + E n _i — 2E n , (12) 

where E n are the total ground state energies. The calculation above finds 
that all odd fillings n = 1, 3, 5 have the same semiclassical pair binding energy 
U 2n+1 = -2S + l-lS-\ 

In Table I we summarize the results for the vibrational contributions to 
the ground state energies and pair binding energies of Cqq u ~ . The semiclas- 
sical results contain the leading three terms in the 1/S expansion. However, 
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one may rightfully worry about higher order corrections since for Ceo the 
experimental estimate is S— 0.2-0.4. In comparing 04Brien's exact results 
for n = 1 [|TTJ to the semiclassical expression for E\ given in Table 1 , we find 
that for S > 0.25, the error in the semiclassical approximation is bounded by 
0.2hu). For S*<0.25, the semiclassical error diverges rapidly. Using S* = 0.25 
for buckminsterfullerene, the semiclassical result is that t/3 is dominated by 
the rotational energy, which enhances the pair binding energy from its classi- 
cal value by the significant amount of hu/3. We must beware however, that 
S = 0.25 is pushing the semiclassical expansion quite far. Since other impor- 
tant interactions have not been considered here (e.g. intermolecular hopping 
and electron-electron interactions[|14|), I refrain from inferring quantitative 
predictions for T c . The results suggest that further investigations of (^) and 
its extensions would be worthwhile. 

In conclusion, I have shown that the dynamical Jahn teller effect in Ceo n ~ 
involves several interesting features. For n = 1,2,4,5 the molecule distorts 
unimodally, giving rise to a pseudo-angular momentum spectrum, plus three 
harmonic oscillators. For n = 3, there is a bimodal distortion, which generates 
a spectrum of a symmetric top rotator, plus two harmonic oscillators. The 
Berry phases of the electronic wave functions have been calculated. They 
determine the allowed pseudo-angular momenta quantum numbers. It would 
be interesting if further spectroscopic investigations of Ceo ions could resolve 
the special structure predicted by the Eqs.([7|,|9]). 
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Table I 



n 


(-2-n, ^n) 


(ni,n 2 ,n 3 ) 


E n /(hu) 


£/ n /(M 



1 
2 


(0,0) 
(9,0) (van) 
(2<7,0) (uni) 


(0,0,0) 
(0,0,1) 
(0,0,2) 


5 

3 2 1 1 
—S +2+6^ 

-45+| 


-25 + 1- 
45-I + 15- 1 


3 
4 
5 
6 


(19,^9) (bi) 
(-20,0) (uni) 
(-9,0) (uni) 
(0,0) 


(1,0,2) 
(2,2,0) 
(1,2,2) 
(2,2,2) 


-35 + 1 + 15- 1 

-45 + 1 
-5 + 2 + g5 

2 


-25 + 1- §5^ 
45-I + 15- 1 
-25 + 1- IS' 1 
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